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ABSTRACT 


The  optimal  detection  of  a  three-dimensional  moving  target  calls  for  the  classical  tech¬ 
nique  of  matched  filtering.  If  a  target  is  modeled  as  a  moving  point  source  with  unknown 
velocity,  then  the  velocity  alone  determines  the  shape  of  the  observed  signal.  Thus,  target 
velocity  is  a  parameter  that  completely  characterizes  the  matched  filter.  We  will  designate 
these  types  of  matched  filters  to  be  the  assumed  velocity  filters  (AVFs)  to  emphasize  the 
velocity  parameter. 

Like  most  matched  filtering  techniques  where  the  signal  parameters  range  in  a  contin¬ 
uum,  the  AVF  must  be  implemented  suboptimaUy  by  quantizing  the  velocity  space.  In  this 
report  we  will  use  a  loss  factor  that  measures  the  average  loss  of  signal-to-noise  ratio  (SNR) 
at  the  output  of  the  matched  filter  due  to  mismatch  of  filter  parameters.  The  loss  factor 
can  be  used  as  a  criterion  for  partitioning  the  velocity  space.  We  will  show’  that,  with  a 
fixed  loss  factor,  the  number  of  filters  required  for  coverage  increases  linearly  as  the  span  of 
the  two-dimensional  velocity  space  increases  quadratically.  The  rate  of  increase  is  further 
reduced  w’hen  the  loss  factor  is  made  proportional  to  expected  target  angular  speed. 


in 


CONTENTS 


Abstract .  iii 

List  of  Illustrations .  vii 

1  INTRODUCTION .  1 

2  THE  DESIGN  OF  THE  ONE-DIMENSIONAL  AVF .  4 

2.1  The  Rectangular-Shaped  Signal .  4 

2.2  The  Gaussian-Shaped  Signal  .  13 

3  THE  DESIGN  OF  THE  TWO-DIMENSIONAL  AVF  .  17 

4  THE  VARIABLE  LOSS  FACTOR .  23 

5  SUMMARY .  29 

Acknowledgment .  29 

References .  30 


v 


LIST  OF  ILLUSTRATIONS 


Figure  Page 

No. 

1  AVF  with  Mismatched  Velocities .  5 

2  Cross  correlation  Function .  8 

3  vp/vp  -  Rectangular  Signal .  11 

4  Number  of  Filters  Required  -  Rectangular  Signal .  12 

5  vp/vp  -  Gaussian  Signal .  15 

6  Number  of  Filters  Required  -  Gaussian  Signal .  16 

7  Filter  3-dB  Loss  Contours  in  Speed-Angle  Space .  21 

8  Fixed  Loss  Factor  (a)  Placement  and  (b)  Growth  Rate  of  Filters .  22 

9  Sensor-Target  Geometry .  24 

10  Variable  Loss  Factor  (a)  Placement  and  (b)  Growth  Rate  of  Filters  ....  28 


Vll 


ON  THE  DESIGN  OF  SUBOPTIMAL  MATCHED 


FILTERS  FOR  THREE-DIMENSIONAL 
MOVING  TARGET  DETECTION 

1  INTRODUCTION 

A  point  target  at  location  x  traveling  at  velocity  v  (see  *)  makes  a  streak  on  the  focal 
plane  of  a  telescope.  The  optimal  detection  of  the  three-dimensional  target  streak  calls  for 
the  classical  technique  of  matched  filtering.1"3  Assume  that  the  focal  plane  consists  of  an 
array  of  solid-state  photodetectors  such  as  the  charge-coupled  devices(CCDs).  Then  we 
can  model  the  image  received  at  the  m-n th  photodetector  cell,  denoted  by  rmn(x,  r),  to  be 
the  sum  of  three  independent  Poisson  processes  —  the  photoelectrons  produced  because 
of  the  target  image,  the  photoelectrons  produced  becaused  of  the  sky  background  noise, 
and  the  electrons  produced  by  the  photodetector  noise  (e.g.,  the  CCD  dark-current)^  The 
rmn(i,u)  itself  is  therefore  a  Poisson  process  (due  to  the  reproductive  property  of  the 
Poisson  distribution)  with  the  average  number  of  electrons  given  by  A,smn(x,  v)  +  A0,  where 
\g  is  the  electron  rate  per  unit  intensity  of  the  target  signal,  Aq  is  the  electron  rate  of  the 
combined  noise,  and  smn  is  the  target  signal. 

Using  the  Poisson  distributions,  the  likelihood  ratio  is  found  to  be 


*  Unless  stated  otherwise,  all  velocities  refer  to  angular  velocities  measured  in  pixels  per  frame  time. 

f  Notably  absent  is  readout  circuit  noise,  which  is  considered  negligible  due  to  recent  engineering 
improvements. 
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A 


jt  (AjSmn  “f  Ao)fmn  CXp[— (AjSmn  +  Ao)] 
Agmn  exp(— A0) 

/  A  \  rmn 

=  nm>n  ^1  +  A™n  j  exp(— A,smn)  . 
Taking  logarithms  of  both  sides,  we  get  the  log  likelihood  ratio 

mn  '  0  /  mn 

Assume  that  AJ<smn  <C  Aq.  Then 


(1) 

(2) 


1  ( ^  .  Aj5mn\  A55m 

,nl1+— J“  — 


(3) 


Substituting  Equation  (3)  into  (2)  yields 

ltt  A  =  T  y  ^mn  —  Ao)«Smn 

A0  mn 


(4) 


Equation  (4)  provides  a  formula  for  the  matched-filter  detection  of  the  stressing  case  of  a 
weak  target  streak.3  For  stronger  targets,  detection  [Equation  (4)]  is  suboptimal,  although 
the  increased  signal  strength  has  more  than  compensated  for  that. 

In  Equation  (4),  the  quantity  A0smn  is  not  data  dependent,  and  is  usually  combined 
with  the  detection  threshold.  In  the  following  work,  however,  we  have  found  it  useful  to 
retain  this  term. 

Let  us  assume  that  the  focal  plane  is  sufficiently  large  and  that  the  detectors  are  suffi¬ 
ciently  small,  and  that  we  can  replace  Equation  (4)  by  an  integral 


In  A  =  ^  /  [r(x)  -  A0]s(x)(fx  .  (5) 

Ao  Jk2 

In  Equation  (5),  x  is  a  two-dimensional  vector  and  R2  is  the  second  Cartesian  product  of 
the  real  line.  We  resort  to  the  integral  form  because  the  mathematics  involved  is  somewhat 
simpler. 
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In  the  case  of  N  statistically  independent  image  frames  with  frame  period  T,  Equa¬ 
tion  (5)  can  be  modified  as 

A  N  r 

In  A  =  Y2  /  (r*(^)  “  ^oM®.  -  ivFT)dx  (6) 

where  we  shift  the  matched  filter  by  an  assumed  velocity  vp  (the  filter  velocity),  correlate 
it  with  the  corresponding  image  frame,  and  then  sum  the  correlations  computed  from  the 
individual  frames. 

Let  uj  be  the  target  velocity.  Taking  the  expected  value  of  Equation  (6)  and  using  the 
mean  of  the  observation 

£[ri(x)]  =  A 8s(x  -  ivrT)  +  A0  (7) 


yields 


namely,  the  autocorrelation  of  the  target  signal.  Obviously  Equation  (8)  is  maximized  when 
Vp  =s  vT.  The  result  has  shown  that  a  maximum  likelihood  estimator  would  choose,  from 
all  possible  velocities,  the  Vp  that  maximizes  Equation  (6),  i.e., 


vT  =  max 
2 If 


(9) 


Notice  that  the  formulation  allows  us  to  match  filter  a  random  signal  by  its  conditional 
mean. 
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2  THE  DESIGN  OF  THE  ONE-DIMENSIONAL  AVF 


We  now  investigate  the  output  SNR  of  the  matched  filters  when  the  assumed  velocity 
vp  deviates  from  the  target  velocity  vT  ( vF  ^  We  start  from  the  one-dimensional  case 
because  it  is  both  simple  and  enlightening.  We  will  divide  the  discussion  of  one-dimensional 
assumed  velocity  filters  (AVFs)  into  two  parts;  in  the  first,  we  assume  a  rectangular-shaped 
target  signal,  and  in  the  second,  a  Gaussian-shaped  target  signal.  The  two  results  share 
similar  characteristics.  Therefore,  owing  to  its  mathematical  tractability,  we  will  apply  the 
Gaussian  approximation  to  the  design  of  the  two-dimensional  AVF  in  the  next  section. 

2.1  THE  RECTANGULAR-SHAPED  SIGNAL 

Suppose  vp  is  mismatched  to  vj  (vp  ^  t>j).  By  interchanging  the  summation  and  the 
integral  in  Equation  (6)  and  by  changing  a  variable,  wre  obtain  the  equivalent  expression 

]nA  =  +  ivfT)  -  A0]|  s(x)dx  .  (10) 

The  equation  suggests  that  we  remove  the  mean  of  the  noise  from  each  image  frame,  stack 
the  frames  according  to  the  filter  velocity  vp ,  and  then  correlate  the  stacked  frame  with 
the  filter  signal.  In  Figure  1  we  have  plotted  the  means  of  three  frames  of  images,  and  we 
have  shown  that  when  vp  is  mismatched  to  t>j,  the  mean  of  the  stacked  signal  will  have  the 
shape  of  a  wedding  cake.  A  loss  of  SNR  ratio  is  thus  incurred  w’hen  the  stacked  signal  is 
correlated  with  the  original  rectangular  signal. 

Let  the  stacked  signal  be  denoted  by  y(x).  Its  mean,  y(x),  can  be  represented  by 
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Figure  1.  AVF  with  mismatched  velocities. 
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(11) 


N  N 

y(x)  =  A5  s(x  -  i*>T T  +  ivFT)  =  £  *(*  +  t'At>T) 

i=0  i=0 

where  Av  =  vF  —  vt- 

Now,  for  convenience  of  mathematical  manipulation  we  assume  that  T  is  sufficiently 
small  and  N  is  sufficiently  large  such  that  the  summation  in  Equation  (11)  can  be  replaced 
by  an  integral.  Then 


r 

y(x)  =  Xa  /  s(x  +  Avt)dt  =  Xs  /  w(t)s(x  +  Avt)dt  (12) 

J-NT/2  JR 

where  R  is  the  real  line  and  w(t)  is  a  rectangular  weighting  function  defined  to  be 

/  t  \  1  ,  \*\  <  NT/2 

«>(<)  =  rcct  = 


(13) 


0  ,  \t\  >  NT/2 


According  to  Equation  (10)  we  will  correlate  y(x)  with  a  filter  signal,  which  we  will  refer 
to  as  sF(x).  We  let  the  target  signal  s(:r)  and  the  filter  signal  s/r(^)  he  rectangular  with 
lengths  vjT  and  vpT,  that  is, 


s(z)  =  rect 


sjr(:c)  =  rect 


(vTt)  ~ 


1 , 

|a*|  <  vpT/2 

0 , 

|a:|  >  vpT / 2 

1 , 

|a:|  <  vFT/2 

0 , 

M  >  vpT/2 

(14) 


(15) 


The  mean  matched-filter  output,  denoted  by  !>mf,  is  given  by 


>m/  = 


=  A 


Aa  /  /  w(t)s(x  +  Avt)sp(x)dx  dt 
JR  JR 

f  w(t)RSiaJAvt)dt 
Jr 

=  zi; /*'“•( HZ&)  x-'W* 


(16) 
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The  third  equality  in  Equation  (16)  is  due  to  the  change  of  variable  x  =  A vt.  RSfSF(x)  is 
the  cross  correlation  of  the  target  and  filter  signals,  pictorially  depicted  in  Figure  2.  The 
rectangular  window  function  w(x/N  AvT)  is  also  plotted  in  Figure  2  in  dashed  lines. 

To  evaluate  Equation  (16),  we  need  to  find  the  area  under  the  trapezoid  and  enclosed 
in  the  window  boundaries,  shown  in  Figure  2  as  a  shaded  area.  We  assume,  without  loss  of 
generality,  that  Av  >  0.  Furthermore,  we  assume  that 


<*  -  1>A’-r  <  vtt 


(17) 


The  inequality  [Equation  (17)]  imposes  the  constraint  that  the  temporal  wundow  is  short 
enough  such  that  the  shifted  target  images  remain  mutually  correlated.  In  Figure  2,  we 
see  that  Equation  (17)  requires  that  the  window  boundaries  be  enclosed  in  the  trapezoid. 
Under  these  assumptions,  it  can  be  showm  that  the  filter  velocity  satisfies 


VT  <  Vf  < 


AT  +  1 
N  -  1 


vt 


(18) 


and  that  the  mean-filter  output  signal  is  equal  to 

(N  -  l)2Ar 


Smf  — 


N  vt  — 


(19) 


Of  equal  interest  is  the  noise  variance  at  the  matched-filter  output.  It  is 


70,m  / 


=  J J J J  u(r))n(x  + vpri,ri)sF{x)w{C)n(y+ VFi,Z)sF(y)dxd.ydridZ 

=  Ao  J  w*(T])dT]  J  s2F(£)dt 

=  X0  NvpT2  .  (20) 


The  ratio  of  the  square  of  the  mean  filter  output  signal  and  the  noise  variance  is  the  SNR 
at  the  output  of  the  matched  filter.  We  denote  it  by  SNRoUt.  Then 
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NAvT  0  AvT  (vT  +  vf)T 


2  2  2 

Figure  2.  Cross  correlation  function. 
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SNR, 


out  — 


C0,mf 

A t[Nyr  -  (N  -  l)2Ar/4]2r2 


XqN  vp 


=  SNR, 


Nvj 

vp 


1  - 


(N  -l)2Av 


4Nvj 


(21) 


where  SNRtn,  the  SNR  prior  to  the  matched  filter,  is  defined  as  the  square  of  the  mean  of 
the  signal  divided  by  the  variance  of  the  noise,  viz . 

A? 


SNR  . 

Dividing  both  sides  of  Equation  (21)  by  SNR,n  yields 


(22) 


SNRguf  _  Nvl 


SNR,-, 


t'F 


1  - 


(N  -  l)2At> 


1  2 


AN  vt 


(23) 


If  the  filter  is  precisely  matched  to  the  signal  (Ar  =  0),  Equation  (23)  is  reduced  to 


SNR, 


out 


SNR, 


=  NvtT2 


Av=0 


(24) 


This  is  the  maximum  gain  in  SNR  that  the  matched  filter  is  able  to  achieve.  Substituting 
Equation  (24)  back  to  (23)  yields 

SNR0U(  SNRout 


SNR,, 


SNR,-, 


Vp 

Av=0  VF 


1  _  (Ar  -  l)2Ar 
4Nvj 


n  2 


(25) 


Equation  (25)  describes  the  average  loss  of  SNR  due  to  velocity  mismatch.  In  the  following, 
we  will  show  that  Equation  (25)  can  be  used  to  design  the  matched  filters  such  that  the 
average  SNR  loss  is  no  greater  than  a  prescribed  loss  factor  p  where 


Vt 

P  ~  vf 


1  - 


(. N  -  l)2Av 
4Nvt 


(26) 
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Let  x  be  vp/vj,  a  be  (N  —  1)2/4JV,  and  /?  be  (N  +  l)2/4 N.  Equation  (26)  can  be 
rewritten  as  a  quadratic  equation  in  x 

«2X2  -  (2a/?  +  p)x  +  P2  =  0  •  (27) 

Equation  (27)  has  a  solution 

X  =  ^  =  2^2  (2q/?  +  P~  ^4q/?P  +  p2)  (28) 
that  satisfies  the  constraint  [see  Equation  (18)].  This  solution  gives  the  largest  vp  that  sat¬ 
isfies  the  given  p,  To  find  the  smallest  vp  that  satisfies  the  same  p,  one  simply  interchanges 
the  roles  of  the  target  and  filter  signals  in  Equation  (26).  The  minimum  ratio  vp/vp  is 
equal  to  1/x,  as  defined  in  Equation  (28). 

The  design  of  the  suboptimal  AVF  amounts  to  using  the  relationship  in  Equation  (28) 
to  partition  the  velocity  space  such  that  the  incurred  average  loss  in  filter  output  SNR  is  no 
larger  than  p.  Two  useful  plots  are  given  in  Figures  3  and  4.  In  Figure  3  we  have  plotted 
Equation  (28)  for  a  set  of  frame  numbers.  For  example,  if  N  —  10  then  the  maximum  vp 
that  satisfies  p  —  0.5  (or  3-dB  loss)  is  vp  =  1.124 vp,  as  indicated  in  Figure  3  by  an  open 
circle.  The  minimum  vp  is  given  by  vp  =  0.89 vp. 

A  useful  property  that  has  surfaced  from  the  derivation  is  that  the  maximum  (as  well  as 
the  minimum)  vp  that  satisfies  a  given  p  is  proportional  to  vp.  That  is,  the  spacing  between 
filter  velocities  increases  sis  the  expected  target  velocity  increases.  Therefore,  a  relatively 
fewer  number  of  filters  is  needed  to  detect  higher  speed  targets  than  that  needed  for  lower 
speed  targets.  This  is  demonstrated  in  Figure  4,  where  we  have  plotted  the  numbers  of 
filters  required  for  three  velocity  ranges  (see  *,  page  1  for  unit  of  velocity):  1  to  50,  1  to 
100,  and  1  to  200.  The  number  of  frames  is  10.  For  example,  at  3-dB  loss,  Figure  4 
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Figure  3.  vp/vj  -  rectangular  signal. 
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Figure  4.  Number  of  filters  required  -  rectangular  signal. 


12 


88409-4 


shows  that  only  22  filters  are  needed  to  cover  the  velocity  range  1  to  200  (indicated  by  an 
open  triangle),  but  20  of  the  22  filters  are  used  to  cover  the  range  1  to  100  (indicated  by  a 
filled  triangle).  The  property  is  particularly  useful  when  we  get  to  the  problem  of  designing 
two-dimensional  AVF. 


2.2  THE  GAUSSIAN-SHAPED  SIGNAL 


In  anticipation  of  the  two-dimensional  case,  we  also  chose  to  analyze  one-dimensional 
AVF  by  using  Gaussian  approximations.  Let 


s(i)  =  \Zrxp  (-£j) 
«,(<)  =  Vf«p(-s 


(29) 

(30) 

(31) 


where 


2  _  ( vtT ) 

(7  y  — 


o\  - 


= 


12 

(t‘FT)2 

12 

(AT)2 

12 


(32) 

(33) 

(34) 


are  the  variances  that  have  been  chosen  to  match  those  of  the  rectangular-shaped  signals. 
The  multiplicative  factors  in  Equations  (29),  (30),  and  (31)  are  chosen  such  that  the  areas 
under  the  Gaussian  signals  match  those  of  the  rectangular  signals. 

From  Equation  (16)  we  have  learned  that  to  compute  the  filter  output  3^/,  we  need  to 
compute  the  cross  correlation  of  s(x)  and  Sjr(x).  For  Gaussian  signals,  this  is  given  by 


D  /  \  l2aTaF 

R8,8F{X)  =  7= — =  =  =r=€XP 


+  °f) 


2(aT  +  aF) 


(35) 
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Therefore,  Sm/  is 


=  hhw{  jbf)R-^x)dJ 


=  A,,/- 


r  yJo%To£+  Av2al 

and  the  noise  variance  at  the  filter  output  is 

36A0 


^0 ,mf  ~  X°J U,2(7?)dr?  J sf(0^  =  — 


Using  Equation  (36)  and  (37),  it  is  easy  to  see  that 

2AGjGfGw 


out 


SNR, 
SNR,„ 


(?T  +  <Tp+  Av2<T2 

SNR0U< 


SNR,-, 


2vp/ vt 


At/=0  1  +  Vp/vl  +  N2(vF  -  VTy/v$ 


(36) 


(37) 


(38) 


Letting  x  =  vp/vp  and 


_  2  vp/vj 

P  ~  1  +  Vp/Vj  +  At2(uf  -  Vj)2/t'j 

we  obtain  the  following  quadratic  equation 


X 


2 


2AT2  +  2/p 
AT2  +  1  X 


+  1  = 


0 


(39) 


(40) 


A  useful  solution  of  Equation  (40)  is  given  by 
vp  1 

X=  v^=  N2  +  1 


A'2  +  -  -  J 

2  (—  —  l)  N2+  (\-l) 

p  V 

\p  J  Vp2  J  m 

(41) 


Similar  to  Figures  3  and  4,  we  have  plotted  Figures  5  and  6  for  the  Gaussian-shaped  signal, 
using  Equation  (41). 

By  examining  the  equations  and  the  figures,  we  conclude  that  the  Gaussian  approxi¬ 
mation  renders  the  same  filter  characteristics  as  the  original  rectangular  signals,  i.e.,  the 
spacing  between  filter  velocities  increases  as  the  expected  target  velocity  increases.  Owing 
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Figure  5.  vp/vx  -  Gaussian  signal 
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Figure  6.  Number  of  filters  required  -  Gaussian  signal. 
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to  its  mathematical  tractability,  we  will  use  the  Gaussian  approximation  in  designing  the 
two-dimensional  AVF. 


3  THE  DESIGN  OF  THE  TWO-DIMENSIONAL  AVF 


Two-dimensional  assumed  velocity  filtering  may  have  mismatches  either  in  the  speed 
or  direction  of  movement.  In  this  section,  we  use  the  Gaussian-shaped  signals  for  the 
investigation  of  SNR  loss  due  to  such  mismatches.  Finally,  we  will  use  the  results  to  partition 
a  two-dimensional  target  velocity  space. 

Without  loss  of  generality,  we  assume  that  the  target  signal,  described  by  length  vjT 
and  width  /i,  is  aligned  with  the  horizontal  axis.  The  filter  signal,  with  length  vpT  and 
width  ft,  intersects  the  target  signal  at  an  angle  6.  In  the  Gaussian  signals,  we  define  a\ 
and  ajp  as  in  Equations  (32)  and  (33).  In  addition,  to  approximate  the  signal  width  we  use 
a  variance 

cl  =  h2 /12  .  (42) 

Now  the  target  and  filter  signals  are  given  by 


s(x)  =  ^exp 
sf(x)  =  ~exp 


(43) 

(44) 


where  x*  =  [xi  X2].  The  covariance  matrix  of  the  target  signal,  Kj ,  is  the  diagonal  matrix 


Kt  = 


<Tj  0 


0  al 


(45) 


The  covariance  matrix  of  the  filter  signal,  Kp,  is  another  diagonal  matrix  rotated  by  a 
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similarity  transformation 


Kf  = 


cos#  sin#  Op  0  cos#  —sin# 

—  sin#  cos#  0  a\  sin#  cos# 

o\  cos2#  +  o\  sin2#  (cl  —  Cp)  sin#  cos# 
(c\  —  Cp)  sin#  cos#  Cp  sin2#  +  a2  cos2# 
The  cross  correlation  of  s(x_)  and  sf(x)  is  given  by 

n  ,  N  72^/1/1x11  A>|  \  1  ,  F,  wl  1 

R •"& = + avi  exp  (Kt + Kf)  £] 


(46) 


(47) 


Using  Equation  (47)  and  the  temporal  weighting  function,  we  can  compute  the  output  signal 
~§m  /  and  the  noise  variance  as 


Smf  = 


_2 

aO,mJ 


A4  J  w(t)Rs>SF(Avt)dt 

_ 72N/121Ar||A>|A< _ 

*y/\KT  +  Kf\[\I°1  +  Av}{Kt  +  KF)~lAv] 

Ao  J™2(n)dri  Js2F(C)di 
432cwXo  [77,  T 

-75T'/|AfI 


(48) 


(49) 


Av  = 


where 

Vp  cos  #  —  vt 
vpsmO 

Using  the  above  equations,  we  can  express  SNRout  as 

1A>]]A'f]A? 


(50) 


SNRout  =  144 
=  144 


fl 

V  7r  | Kt  +  KF\^\K7\aw[l/al  +  A v<(KT  +  KF)~' Av]A0 

[2  \KT\\KF\  SNR,„ 

V*i 


(51) 


|  Kt  +  KF\^\K^\aw[l/al  +  Av'(KT  +  KF)~lAv ]  * 

Therefore,  we  can  write  the  following  equation  for  the  description  of  the  average  loss  of 
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SNR  due  to  the  velocity  mismatch 


SNR0Ui  SNRout 


SNR,„  SNRi„ 


4x/i rap 


LF I 


Av=o  \Kt  +  I<f WUUol  +  Avt(Kj  +  A'f)-1  Au] 


Now  let 


P=TT7 


4\/fA'j||A> 


|Aj  +  AVK  [l/o-a  +  AuHA't  +  Kf)~'Av]  ’ 


It  can  be  verified  that 


\/\Kt\\Kf\  =  oToFa2h 

\Kt  K f\  =  4“  +  Gj'G'jr  sin20  -}-  CTjCT2  cos2$ 

+  (7p(7^  cos2#  +  <7  4  sin2# 

1 


(AV  +  AV)"1  = 


(52) 


(53) 


(54) 


(55) 


(56) 


(gj  +  <yjr)0fc(  1  +  COS2  #)  -f  (cTjCT^  +  <74)  sin2# 

<7^  sin2#  +  <72(1  +  cos2#)  (crp  -  cr2)  sin#  cos# 

(cT^  -  cr^)  sin #  cos #  Gj  +  g\  cos2#  +  cr2  sin2# 

Using  Equations  (50),  (54),  (55),  and  (56),  and  letting  %  =  vp/vj  and  a  =  h2/vjT 2 
Equation  (53)  can  be  rewritten  as  a  quartic  equation  in  x 


(4pAr2sin2#cos2#)x4  -  (4pAr2sin2#cos#)x3 

+[pa(l  +  cos2#)  +  psin2#  +  2pN 2  sin2#  +  pAr2a(l  +  cos2)  cos2#  -  2pN2a  sin2#  cos2  # 

—2pN2a  sin2#  cos 2 6  +  pN2a  sin4#]x2 

+[-4 a  -  2pAr2a(l  +  cos2#)  cos#  +  2pN2a  sin2#  cos#]x 

+[pa(l  +  cos2#)  +  pa2  sin2#  +  pN2a(l  +  cos2#)] 

=  0  .  (57) 


19 


Equation  (57)  can  be  solved  for  a  given  set  of  6  (see  footnote).  The  results  may  be 
used  to  partition  the  v  -  0  plane.  Keeping  N  =  10,  we  plot  the  filter  3-dB  loss  contours  in 
Figure  7  for  a  set  of  carefully  spaced  filters.  A  target  with  speed  and  angle  contained  in  any 
of  the  cells  will  be  detected  by  the  corresponding  filter  with  an  SNR  loss  <  3  dB.  Because 
of  the  irregular  shape  of  the  cells,  proper  overlapping  wfil 1  be  required  to  fill  the  gaps.  For 
clarity,  we  do  not  show  such  overlapping  in  this  and  all  subsequent  figures. 

The  design  of  the  suboptimal  two-dimensional  AVF  amounts  to  the  partitioning  of  the 
velocity  space  using  Equation  (57).  As  an  example,  we  have  designed  filters  with  a  maximum 
3-dB  loss  for  the  coverage  of  the  tq  -  v2  plane,  where  1  <  yjv\  +  v\  <  10.  The  partitioning 
of  the  velocity  space  is  shown  in  Figure  8(a).  A  very  useful  property  of  the  design  can 
be  noticed  in  the  figure,  that  is,  the  radial  spreading  of  the  resolution  cells.  Had  the  size 
of  the  cells  remained  unchanged,  the  number  of  filters  required  for  coverage  would  have 
increased  quadratically  as  tq  and  v2  increased.  However,  due  to  the  radial  spreading  of  the 
cell  size,  we  expect  this  growth  to  be  much  slower.  A  count  of  the  number  of  filters  (less 
those  overlapping),  shows  a  linear  growth  rate,  as  indicated  in  Figure  8(b).  This  property 
is  crucial  for  making  efficient  implementation  of  AVFs  feasible. 


Two  situations  may  occur  in  solving  Equation  (57): 

1.  There  are  no  real  solutions.  For  the  given  0,  the  loss  factor  p  cannot  be  satisfied. 

2.  There  is  one  and  only  one  pair  of  real  solutions.  The  real  solutions  are  the  solutions  of  interest. 
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Figure  7.  Filter  3-dB  contours  in  speed-angle  space. 
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Figure  8.  Fixed  loss  factor,  (a)  placement  and  (b)  growth  rate  of  filters. 
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4  THE  VARIABLE  LOSS  FACTOR 


In  the  begining  of  this  section,  we  take  a  slight  change  of  course  and  investigate  how 
the  filter-output  SNRotif,  varies  as  a  function  of  target  speed.  Later,  we  will  return  to  the 
topic  of  AVF  design  and  we  will  use  the  result  developed  to  define  a  variable  loss  factor, 
as  opposed  to  the  previously  used  fixed  loss  factor.  Through  the  same  example  of  the  last 
section,  we  will  demonstrate  that  even  greater  efficiency  can  be  achieved  with  the  variable 
loss  factor. 

Assuming  a  perfectly  matched  filter,  Equation  (51)  can  be  simplified  to 


oi™  I  _  *  /6  \]ohawvTT 
SM^out|^v=o  —  6y  jj.  ^ 


(58) 


The  SNRotit  is  seen  to  be  proportional  to  vj,  a  well-known  characteristic  of  matched  filtering. 
Hidden  in  Equation  (58)  is  that  SNRotit  may  also  be  related  to  vj  through  its  association 
with  A,,  the  electron  rate  due  to  the  target  signal.  The  A,  is  inversely  proportional  to  both 
the  target  angular  speed  and  the  square  of  the  distance  between  the  sensor  and  the  target. 
Let  r  be  this  distance,  then 


A,  = 


VjT * 


(59) 


where  C,  the  constant  of  proportionality,  accounts  for  the  visual  magnitude  of  the  sun, 
effective  cross  section  of  the  target,  frame  integration  time,  telescope  aperture,  CCD  photon 
efficiency,  etc. 

Based  on  Equation  (59)  we  will  treat  three  different  cases  of  ground-based  sensors: 
Case  1.  The  sensor  is  located  at  the  center  of  the  earth  [as  shown  in  Figure  9(a)].  By 
Kepler’s  third  law 


6  =  vtT  <x  r~3'2 


(60) 
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Figure  9.  Sensor-target  geometry. 
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Substituting  Equation  (60)  into  (59)  yields 


A,  a  Vj3  .  (61) 

Finally,  substituting  Equation  (61)  back  into  (58)  yields 

SNR ^  oc  4/3  .  (62) 

Case  2.  The  sensor  is  located  on  the  surface  of  the  earth  and  the  target  is  directly 
overhead  [as  shown  in  Figure  9(b)].  For  simplicity,  we  ignore  the  rotation  of  the 
earth.  From  the  figure,  Kepler’s  third  law  is  stated  as 

0  a  (r  +  te)"3/2  (63) 

where  V£  is  the  local  radius  of  the  earth.  Since 

rp  fs  d  (r  +  r^)-1  (64) 

the  linear  velocity  d  relates  to  r  through 

d  oc  (r  +  te)-1'2  (65) 

and  then  the  angular  velocity  vj  relates  to  r  through 

6  =  vtT  oc  r~l(r  +  ve)~1/2  cost) 

~  r-1(r  +  rE)~1/2  (66) 

where  we  have  assumed  that  the  angle  8  is  small.  We  consider  the  following  two 
extremes: 

I.  When  r  >  r^,  then 

vt  oc  r_3/2  .  (67) 
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Since  Equation  (67)  is  identical  to  (60)  in  Case  1,  we  conclude  that  Equa¬ 
tion  (62)  holds  in  this  case  too. 

II.  When  r  <C  te ,  then 

vt  oc  r”1  .  (68) 

Substituting  Equation  (68)  into  (59)  yields 

A*  oc  vt  .  (69) 

That  is,  Xs  increases  as  the  angular  target  speed  increases.  Therefore,  in 
this  case 

SNRovt  vt  •  (70) 

Case  3.  The  sensor  is  located  on  the  surface  of  the  earth.  The  target,  sensor,  and  center 
of  the  earth  form  an  angle  (3  [as  shown  in  Figure  9(c)].  This  is  a  generalization 
of  Case  2.  Applying  Kepler’s  third  law  again  we  have 

rp  oc  iT3/2  (71) 

d<xR~1'2  (72) 

and 

6  =  vjT  a  r~1R~1^2  cos  a  .  (73) 

By  the  law  of  cosines  we  can  find  cos  a.  It  is  given  by 

cos  a  =  (r  —  r£Cos/?)i?-1  .  (74) 

Substituting  Equation  (74)  into  (73)  yields 

6  =  vjT  oc  r-1(r2  +  r%  —  2rr£Cos/3)_3/,4(r  —  re  cos/?)  .  (75) 
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Now  it  can  be  verified  easily  that  the  two  extremes  of  Case  2  hold  approximately 
true  even  when  the  target  is  not  directly  overhead. 

In  all  three  cases  of  ground-based  sensors,  we  have  shown  that  SNRout  is  at  least  pro¬ 
portional  to  v^/3.  The  implication  of  the  above  relationships  is  that  fast  targets  are  brighter 
and  therefore  they  are  easier  to  detect  (greater  probability  of  detection  at  the  same  prob¬ 
ability  of  false  alarm)  than  slow  targets.  For  the  purpose  of  efficient  utilization  of  detector 
resources,  it  may  be  desirable  to  design  a  detector  with  uniform  sensitivity  to  various  target 
speeds.  This  can  be  achieved  through  adjusting  the  loss  factor  in  anticipation  of  different 
expected  target  angular  speed. 

To  demonstrate,  we  will  adopt  a  conservative  linear  relationship.  Let  Vq  be  the  lowest 
angular  speed  the  system  is  designed  to  detect;  let  v  >  t>o  be  that  of  an  arbitrary  target; 
and  let 

.  (76) 

Letting  v0  =  1,  we  rework  the  example  in  the  last  section.  The  loss  factor  pf  is  now 
l/2v  for  p  =  1/2.  The  partitioning  of  the  velocity  space,  based  on  the  solutions  of  Equa¬ 
tion  (57),  is  plotted  in  Figure  10(a).  Notice  that  the  size  of  the  resolution  cells  increases 
quite  dramatically  compared  with  that  of  Figure  8(a).  Figure  10(b)  illustrates  the  much 
slower  growth  rate  demonstrated  in  this  example.  For  comparison,  38  filters  are  needed  to 
cover  the  10  pixels  per  frame  velocity  space  (with  no  overlapping  of  cells)  in  this  example, 
whereas  191  filters  were  needed  when  the  loss  factor  was  held  constant. 
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Figure  10.  Variable  loss  factor,  (a)  placement  and  (b)  growth  rate  of  filters . 
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5  SUMMARY 


The  problem  of  matched-filter  detection  of  three-dimensional  moving  targets  is  consid¬ 
ered.  Both  the  signal  and  noise  are  assumed  to  be  Poisson  processes.  The  matched  filters  are 
derived  based  on  a  weak  signal  model.  They  are  called  the  assumed- velocity  filters  (AVFs) 
because  velocity  is  the  only  filter  parameter.  For  efficient  implementation  of  the  AVF  it  is 
necessary  to  partition  the  velocity  space.  For  this  purpose,  we  have  introduced  a  loss  factor 
which  characterizes  the  loss  of  SNR  due  to  velocity  mismatches.  For  a  prescribed  maximum 
loss,  we  have  derived  mathematical  equations  that  can  be  used  to  compute  the  boundary 
of  velocity  resolution  cells.  The  cells  can  then  be  used  to  partition  velocity  space.  We  have 
shown,  through  simulation,  that  the  number  of  matched  filters  required  increases  linearly 
as  the  span  of  the  velocity  space  increases  —  a  desirable  property.  Finally,  we  discussed  the 
use  of  the  variable  loss  factor  where  wTe  have  explored  the  increase  in  signal  intensity  due  to 
the  closeness  of  a  target.  Using  the  variable  loss  factor,  the  number  of  filters  required  for 
coverage  is  further  reduced  to  only  a  few  dozen  for  single  quadrant  coverage. 
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